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Abstract. A detailed knowledge of the Galactic radio continuum is of high interest for studies of the dynamics 
and structure of the Galaxy as well as for the problem of foreground removal in Cosmic Microwave Background 
measurements. In this work we present a full-sky study of the diffuse Galactic emission at frequencies of few GHz, 
where synchrotron radiation is by far the dominant component. We perform a detailed combined analysis of the 
extended surveys at 408, 1420 and 2326 MHz (by Haslam et al. 1982, Reich 1982, Reich & Reich, 1986 and Jonas 
et al. 1998, respectively). Using the technique applied by Schlegel et al. (1998) to the IRAS data, we produce 
destriped versions of the three maps. This allows us to construct a nearly-full-sky map of the spectral index and 
of the normalization factor with sub-degree angular resolution. The resulting distribution of the spectral indices 
has an average of j3 = 2.695 and dispersion ap = 0.120. This is representative for the Galactic diffuse synchrotron 
emission, with only minor effects from free-free emission and point sources. 

Key words. Galaxy: general - Methods: data analysis. 



1. Introduction 

The detailed study of radio continuum emission allows 
a direct evaluation of important parameters describing 
the dynamics and structure of the Galaxy, such as mag- 
netic field distribution and electron energy distribution. 
Moreover, diffuse Galactic radio emission is one of the 
main sources of unwanted signal in observations of the 
Cosmic Microwave Background (CMB) at v ;$ 50 GHz. 
High precision imaging of the CMB fluctuations is a 
unique cosmological tool to determine fundamental pa- 
rameters that describe the structure and evolution of the 
Universe and directly probes the early phase of cosmic 
history. 

The new generation of CMB imaging experiments (in 
particular the WMAP 1 and PLANCK 2 space mis- 
sions) call for extended, detailed information of Galactic 
foreground emissions. While such exquisite precision mea- 
surements provide a great information on foreground emis- 
sions at microwave frequencies, a deep knowledge at radio 
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and sub-millimetre frequencies of all emissions that are 
between us and the Last Scattering Surface where CMB 
photons originate significantly contribute to the accurate 
exploitation of microwave data. In particular great atten- 
tion has to be devoted to Galactic foreground emissions 
with their own power spectra and frequency dependence. 

The three main physical processes that produce 
Galactic diffuse emissions are synchrotron, free-free (or 
thermal bremsstrahlung) and dust emission. While free- 
free does not dominate in any observed frequency band, 
synchrotron and dust are the major contributors at low 
and high frequencies, respectively. Therefore, it is conve- 
nient to deal with these two foreground emissions using 
template maps obtained from observations where such 
emissions overwhelm the others. In this way it is possi- 
ble to cross-correlate experimental CMB data with exter- 
nal data sets ( |Banday et al. 20031 Bennett et al. 2003) to 
clean CMB signal from foreground signal as well as per- 
forming component separation with the use of priors on 
foreground emissions UMaino et al. 2003|l . 

As for the dust emission is concerned, a detailed work 
has been done by Schlegel, Finkbeiner & Davis (1998, 
hereafter SFD) to combined IRAS with COBE-DIRBE 
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data in order to obtain a template for thermal dust emis- 
sion with IRAS angular resolution and COBE-DIRBE 
calibration quality. 

In the case of synchrotron emission, observations suffer 
from many systematic effects which depend on the spe- 
cific experimental setup and observing strategies. Due to 
the need of covering the whole sky, many instruments per- 
form observations that will be finally combined in a unique 
data set. This introduces problems such as non-uniform 
zero-levels, resulting in a normalization of the signal when 
combining different patches of observation taken by differ- 
ent instruments located in different places. Furthermore 
the involved technology makes use of radiometer detec- 
tors which typically show 1/ / noise induced by gain fluc- 
tuations that leads to striping effects in the maps when 
combined with the observing strategy. 

In this paper we address the issue of cleaning low fre- 
quency maps from their striping effect employing a tech- 
nique similar to that used by SFD on the IRAS plates. 
We apply this to the three largest low frequency maps: 
the 408 MHz map by Haslam et al. (1982), the 1420 MHz 
map from Reich (1982) and Reich & Reich (1986), and the 
2326 MHz map by Jonas et al. (1998). 

Our goal is to produce reliable templates of spectral in- 
dex and normalization factor of the power law description 
of the Galactic emission at « GHz frequencies, by remov- 
ing the main observation artifacts and taking into account 
the survey calibration and zero level uncertainties. As the 
recent WMAP results show, any template of the Galactic 
emission is frequency dependent l|Bennett et al. 2003)1 . A 
comparison between the synchrotron spectral index map 
at WMAP frequencies and that derived from the analysis 
of radio surveys at « GHz frequencies will add informa- 
tion on the physical processes of electron energy loss and 
the properties of Galactic structures. 

This paper is organized as follows. In Sect.Elwe briefly 
summarize the main properties of the considered surveys. 
The description of the first steps of the adopted data anal- 
ysis method, source identification and destriping, is given 
in Sect. |21 while in Sect. 0] we discuss the minimization 
of the discontinuities of the maps of spectral index and 
normalization factor of the radio emission power law due 
to the survey calibration and zero level errors and present 
our results. Finally, we discuss the results and draw our 
main conclusions in Sect. [5] 

2. Data description 

Synchrotron radiation, arising from cosmic ray electrons 
spiralizing into the interstellar magnetic field, is dominant 
over bremsstrahlung (free-free) emission up to 10-20 GHz. 
The emission spectrum has a power law shape, whose spec- 
tral index varies across the sky reflecting variations of the 
electrons density and magnetic field; a realistic model of 
the synchrotron emission has to take into account this 
significant spatial changing, besides the steepening of the 
spectrum with increasing frequency, due to the increase 
of energy losses by high energy electrons. Free-free emis- 



sion comes from electron-ion scattering with a power law 
shape in frequency; since the diffuse emission at rj GHz 
frequency is dominated by synchrotron radiation, no free- 
free templates are available at this frequencies. 

The only full sky map at decimeter range is the well 
known 408 MHz map ( Has lam et al. 198~2|l . Obtained from 
four different experiments, the overall map has an angular 
resolution of 0.85° and a full beam brightness sensitivity 
of 2 K. The estimate of calibration and zero level errors 
quoted by the authors are 10% and 3 K respectively. 
The 1420 MHz map l|Reich 19821 IReich fc Reich 1986|l . 
covering the sky at —19° < 5 < 90°, has a quoted full 
beam brightness sensitivity of 50 mK (0.59° angular reso- 
lution) . The authors estimate 5% of calibration error and 
0.5 K error of the zero level. 

The higher frequency map (2326 MHz. l.Ionas et al. 1998)1 . 
covering the sky region -83° < S < 13° for Oh < RA < 
12h and -83° < S < 32° for 12ft, < RA < 24ft, has an 
angular resolution of 19' and a full beam brightness sen- 
sitivity of 25 mK. Calibration and zero level errors are 
5% and 80 mK, respectively. The receiver was sensitive to 
linearly polarized radiation. 

The combination of gain drifts and observational strat- 
egy results in artificial structures in the three final maps: 
these appear like stripes particularly intense in the scan 
directions. A careful analysis of the systematic effects in 
the 408 and 1420 MHz maps has been performed by Davies 
et al. (1996). Following their considerations, the most evi- 
dent artifacts in the 408 MHz map are the nearly vertical 
stripes at constant right ascension generated by the ele- 
vation scanning strategy and not completely eliminated 
by the nodding scan technique used in the experiment. 
Horizontal stripes are also evident at declination < 20°. 
In the case of the 1420 MHz map the inclined striations 
reflect the azimuthal scanning strategy. Further artifacts 
show up in the horizontal direction mostly in the declina- 
tion range between 20° to 50°. After point source extrac- 
tion, the authors remove the major Fourier components 
which produce horizontal and vertical striations on an- 
gular scales < 5° in the Haslam map and diagonal and 
horizontal striations on scales < 5° and < 9° respectively 
in the 1420 MHz map. As they point out, this technique 
removes the easily visible angular frequency leaving spuri- 
ous structures on scales of several tens of degrees. Vertical 
striations parallel to the scanning direction are visible in 
the 2326 MHz map too and remain even after the smooth- 
ing applied by Jonas et al. (1998). 

3. Description of the destriping method 

To clean the data from the artifacts mentioned above, it 
is necessary to select the contaminated directions and re- 
move the excess signal; for this purpose we adapted the 
method proposed by SFD to destripe the IRAS map. We 
remember that, although a deep comprehension of the ex- 
perimental/observational origin of stripes in the consid- 
ered survey certainly helps the data exploitation, this and 
other kinds of destriping algorithms (see e.g. Keihanen 
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Fig. 1. The block-diagram summarizes the steps of the destriping procedure. 

Fig. 2. A patch of the 408 MHz map is shown (left) together with the same destriped (right). The lower is the stripe 
image; note the diagonal features together with the dominant vertical striations. For a better qualitative comprehension, 
different color tables have been used for the two upper images and the lower one. For the same reason a different color 
table has been used for the visualization of the entire map (Fig. 0J . 



al. 2003, and references therein, for a different destriping 
approach) have the remarkable property to work almost 
independently from the details of the stripe pattern. For 
a comprehensive explanation of the algorithm concept we 
refer to their original paper; here we explain briefly the 
main steps of the process, that are also summarized in 
the block diagram of Fig. ^ The analysis is not performed 
on the map as a whole, since the code works on "squared" 
patches. We divided the maps of the surveys at 408, 1420, 
and 2326 MHz respectively into 44, 33, and 22 partially 
overlapping patches of 152 x 152 pixels with angular size 
A9 x A<j) = 0.35° x 0.35°. The destriping code acts on 
maps gridded in standard polar coordinates 9, <j) by using 
an Equi- Cylindrical Projection (ECP) pixelization, i.e. a 
rectangular, not equal area, projection of the sphere. 
The algorithm first performs a point source removal, se- 
lecting regions whose signal exceeds the median filtered 
image by a given threshold. The source signal is then re- 
placed with the median signal from an anulus around the 
source (SFD); this allows destriping to be performed in 
the Fourier space, where power from bright sources may 
add and overwhelm many wavenumbers. 
The Fourier domain, parameterized in polar coordinates 
(k r ,kg), is then divided into 90 bins of 2° in kg and a back- 
ground power Pback is calculated by "finding the median 
at each 6 and then median filtering the result" (SFD). The 
contaminated wavenumbers are those for which the power 
ratio 70 between the power, Pg, in each kg bin and the 
background power, Pback, differs significantly from unity. 
The difference, APg, between the power as a function of 
9 and a mean "good" power, P avc (obtained by averag- 
ing all power at each (K r , Kg) with jg lower than the 
stripes threshold 7), is then inversely transformed to get 
the image of stripes only, to be subtracted from the orig- 
inal patch. 

The original resolution of the 2.326 GHz map is approxi- 
mately the same as the pixel size used in this analysis; we 
tested the influence of a better sampling on the destriping 
procedure adapting the algorithm to the same map with 
2048 x 1024 pixels (each pixel of 0.175° x 0.175°). The re- 
sults show that the pixel by pixel temperature difference 
between the same destriped patches with different pixel 
size is, on average, less than 5% for 70-80% of the pixels. 
Given that the three frequencies maps will be jointly used 
in the following analysis, we choosed to perform the de- 
striping on the 2.326 GHz map with the same pixel size 
as the others. 

The real space image of a typical Haslam patch is shown 
in Fig. together with the same patch after destriping. 



In the same figure the real space image of the subtracted 
stripes is shown. The destriped patches are then joined 
together to obtain the destriped map 3 . 

3.1. Determination of the thresholds 

Criteria for the choice of the source extraction threshold 
h s and the stripes threshold 7 have to be carefully ad- 
dressed in order to optimize the destriping procedure. We 
performed several tests to determine the best thresholds 
for source and stripe removal and to estimate the relative 
errors and the effect on the destriped maps. The choice 
of the thresholds does not have the same effect on every 
patch of one map. However, the values of h s and 7 have 
been mantained fixed for all the patches of the same map, 
since the division of each map in patches is useful in the 
analysis but does not reflect a real experimental situation. 
We therefore look for the values that give the most uni- 
form results on the whole map. 

A first test for the effectiveness of the source removal 
procedure is to look at the cleaned image to check for the 
presence of residual sources. It should be noted that the 
optimal choice for stripe threshold 7 is strictly related to 
source removal. Reminding that the goal is to obtain an 
image of stripes to be subtracted to the original patch, 
two opposite non-optimal choices of the source threshold 
h s could affect the destriping method. The choice of a 
high threshold h s (at limit, avoid source subtraction at 
all) implies that power from unsubtracted sources could be 
attributed to stripes and then uncorrectly removed from 
the original map. On the other hand, the choice of a too 
low threshold h s implies that not only the signal from 
point sources is subtracted but also power from stripes 
that will not be removed in the destriping procedure. 

To optimize the stripe removal, in the Fourier space 
we first eliminate the peaks higher than 7 as a function 
of the direction Kg (see Fig. |3J). Stripes have a high level 
of directionality and this gives a good criterion for the 
choice of the threshold value: 7 has to be not too low in 
order to maintain this feature and avoid subtraction of 
real astrophysical signals. On the other side, a too high 7 
would eliminate only stripes with very high signal, often 
in one single direction. Fig. |31 shows the pattern of the 
ratio Pg I Pback as a function of the angle 9 for the same 
patch of Fig. 



3 For three patches of the Haslam map the procedure failed 
in finding a lower signal anulus to replace the source signal; in 
these cases the algorithm did not perform destriping. 
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Fig. 3. Ratio between the power, Pg, in each kg bin and the background power, Pback, for the same patch of Fig. EI 
as a function of the angle 8. Different thresholds 7 are shown (horizontal dashed lines). The highest value selects only 
very intense stripes; by lowering the threshold, lower peaks are also included. 





ST < 1% 


1% < ST < 5% 


5% < ST < 10% 


ST > 10% 


408 


36.4% 


55% 


7.0% 


0.6% 


1420 


89.14% 


10.7% 


0.15% 


0.01% 


2326 


18.3% 


50.6% 


22.5% 


8.6% 



Table 1. Distribution of percentage temperature variation ST values in four different ranges and the corresponding 
pixel percentage. ST is evaluated as the difference after and before destriping. 



After the identification of the optimal stripe thresh- 
old 7, we tested the influence of this choice on the de- 
striped map, comparing it with the same map destriped 
with a higher value of 7. By changing 7 by 20% (i.e. to 
a value well beyond a reasonable choice of the threshold) , 
the number of pixels for which the temperature changes 
more than 10% are less than ~ 1% in the 408 MHz map, 
~ 0.1% at 1420 MHz, and reaches a few percent in the 
2326 MHz map (each fraction is referred to the sky cov- 
erage of each survey). As shown in Tabled the destriping 
procedure does not change the intensity level of the sur- 
veys except for a small percentage of the pixels for the 408 
and 1420 MHz maps and for less than 10% of the pixels 
for the 2326 MHz map. 

3.2. Resulting maps 

The resulting maps are shown in Figs. 01 03 HO We also 
produced the three destriped maps without point sources 
using the same technique that has been implemented in 
the destriping algorithm. The technique works well for iso- 
lated sources; in the case of extended structures, however, 
the replacement of the high source signal with the median 
background cause the smoothing of the structure borders. 
We thus preferred to leave the sources in the maps; their 
contribution to the spectral index distribution will be dis- 
cussed in Sect. 14.31 Maps are also available in HEALPix 4 
format (Gorski et al. 1999) with A^idc = 256. The conver- 
sion from ECP to HEALPix maps is performed working 
in real space via an intermediate step, i.e. an HEALPix 
map with higher resolution (N s id c = 4096, the maximum 
available within HEALPix packages). For each HEALPix 
pixel of the higher resolution map, the conversion code 
finds the corresponding ECP pixel and assign its signal 
value. The map is finally degraded to A^idc = 256. 

4. Spectral index and normalization factor maps 

One of the main goals of this work is the study of the 
spectral index of the Galactic diffuse emission over the 
whole sky, except the region —90° < 5 < —83° where only 
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5 Details on the conversion procedure can be found in 
Cappellini et al. (2003a, 2003b). 



408 MHz data are available. Before evaluating the distri- 
bution of the spectral index /3, we subtracted the unre- 
solved extragalactic background and CMB absolute tem- 
perature at 408 and 1420 MHz following Reich & Reich 
(1988). The unresolved extragalactic background for the 
2326 MHz map has been estimated using the method pro- 
posed by Lawson et al 1987. In addition a baseline offset 
of —0.13 K for the 1420 MHz map has been considered 
(see Table |5J). The two maps at higher frequencies have 
been convolved to match the Haslam map resolution. 

We carefully considered the statistical and systematic 
errors of each surveys in the analysis; in the definition of 
the overall temperature error for each pixel, we include 
the quoted sensitivity, er SC ns, the errors on the zero level, 
(T01, and calibration, er ca i, and a correction factor, c a , to 
account for the different area of the pixels 



c a = L^sm 6»j, 



(1) 



where L is the pixel side in steradian and 61 = n — Lj is 
the colatitude of the position of j — th pixel lower border. 
The convolution of the 1420 and 2326 MHz maps to the 
resolution of the 408 MHz map has an impact on the sen- 
sitivity; we took it into account in the error budget. The 
experiment that produced the 2326 MHz map was sensi- 
tive to a linear polarization direction l|.Tonas et al. 19 98). 
while the other two surveys give the total (non polar- 
ized) intensity. The available polarization surveys at low 
frequency do not cover large regions of those observed 
by the Rhodes/HartRAO 2326 MHz telescope (for ex- 
ample IDuncan et al. 1 997 and 1999) or are undersam- 
pled ( TBrouw fc Spoelstra 1976| ); thus a pixel by pixel 
subtraction of the linear polarized component of the 
2326 MHz map based on the polarization data is not fea- 
sible. We prefer to account for this discrepancy between 
the maps including the conservative term <r po i,2326 = 20% 
in the error evaluation of the 2326 MHz temperature 
(Bacciga lupi et al. 2001| ). The final expressions for the 
overall errors at the three frequencies are: 



(Csens, 408^408 



1/2 



dT- 



1420 



^01,408 + (^4085cal,408) 2 
(CTsons,1420#142o(^142o/^408)/ Ca ) 2 
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Fig. 4. Top: original 408 MHz map HHaslam et al. 1982|l . Bottom: destriped 408 MHz map (this work). 

Fig. 5. Top: original 1420 MHz map HReich 19821 IReich fc Reich 1986jl . Bottom: destriped 1420 MHz map (this work). 



Fig. 6. Top: original 2326 MHz map \ Jonas et al. 1998}) . Bottom: destriped 2326 MHz map (this work). For an easier 
division in relatively large squared patches, the region 13° < 6 < 32° of the 2326 MHz map has not been included in 
the data analysis steps described in Sects. 3.1 and 3.2. 



^01,1420 + (^1420<7cal,142o)' 



1/2 



dT, 



2326 



(CT S cns,2326#2326(#2326/ #408)/ C a) 2 + 



°01,2326 



(T23265cal,2326) 2 + 
2 -| 1/2 

(T'23260'pol,2326) 



(2) 



where h is the full width half maximum of the experi- 
ment at each frequency. Reminding that the relation be- 
tween synchrotron radiation temperature and frequency is 
a power law, 



T v = av 
and 

\ogT v = a - (3\ogv . 



(3) 



(4) 



where loga = a, the spectral index and the normalization 
factor and their errors are simply given by: 



\og{T Vl /T v , 



log(i/i/i/ 2 ) 

log^ilogT^ - log^logT^ 
log^i - log^ 2 



(5) 
(6) 



a -p = logi^i - \ogv 2 \ 



G a = \logV! - logV 2 \ 1 



dT x 



.TiloglO 



T 2 logl0/ 



\ 2 

) (7) 




dTi 



TiloglO 



dT, 



T 2 logl0 



(8) 



where 1 and 2 are two generic frequencies. 

In the central region where the three maps overlap, we 
perform a least square fit to evaluate the (3, a and their 
errors l|Press et al. 19921 assuming a constant (3 between 
408 and 2326 MHz. We repeated the same analysis in- 
cluding only the statistical errors in the error budget. We 
estimated the influence of the stripes threshold choice on 
the spectral index distribution: the mean error is ~ 1%, 
thus negligible in the error budget. 

As we expect, the average spectral index and disper- 
sion are not significantly affected by destriping, while sub- 
stantial changes are produced on a pixel by pixel basis 
(up to A/3 ~ 0.4). Fig.[7|shows the difference A/3 between 
the spectral indices evaluated with destriped and original 
maps. 



4.1. Frequency variation analysis 

The steepening of the synchrotron spectrum with increas- 
ing frequency has to be taken into account in any esti- 
mate of the Galactic radio emission. We tested the as- 
sumption that the spectral index between 408 and 2326 
MHz does not change significantly when evaluated with 
different combinations of the three surveys. In the overlap- 
ping region, where data at the three frequencies are avail- 
able, we compare /?408,i420 with /3i420,2326- Fig. |H1 shows 
the distribution of the difference between the spectral in- 
dex evaluated with 408 and 1420 MHz data, /3408,i420, and 
the spectral index evaluated with the 1420 and 2326 MHz 
data, ,$1420,2326 • The shape of the distribution indicates 
that there is not a clear steepening of (3 when using these 
surveys. Moreover, we compare pixel by pixel the distri- 
bution of the difference with the standard deviation crp, 



where, for each pixel o ^ = a 



'/3(408,1420) t <J ,3(1420,2326) 



only 



for the 5% of the pixels the difference between the spec- 
tral indices is larger than the estimated error. The same 
analysis performed with the statistical errors only shows 
that for the 25% of the pixels the difference is larger than 
C/3,stat- We conclude that, generally, the estimated error 
is larger than the frequency variation of the spectral in- 
dex (3 between 1420 and 2326 MHz. This motivates the 
idea to jointly use, where possible, all the three surveys to 
derive a single map of spectral indices and normalization 
factors in the frequency range 0.408-2326 MHz. As dis- 
cussed in Sect. |3 these surveys are affected by zero level 
and calibration uncertainties; the latter are particularly 
relevant at 2326 MHz. In the next subsection we discuss 
an attempt to reduce, at least in part, the impact of these 
uncertainties in the joint analysis of the three surveys. 

4.2. Discontinuities 

As one might expect, the (3 and a maps show two hori- 
zontal lines of discontinuity corresponding to the borders 
of the maps at 1420 and 2326 MHz. Of course these dis- 
continuities arise from measurements and inter-calibration 
errors and do not have any astrophysical meaning. 
A smoothing of the data to a lower resolution can be eas- 
ily applied to eliminate the appearence of this effect but 
would not eliminate its causes. The same effect has been 
observed by Giardino et al. (2002) who combined the 408, 
1420 and 2326 MHz maps to obtain a large angular res- 
olution (10°) map of j3. To eliminate the discontinuities, 
they tried different techniques looking for a baseline offset 
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Survey 


calibration 


Zero level 


Baseline 


Extragalactic 


Frequency (MHz) 


error (per cent) 


error (K) 


correction 


background (K) (K) 


408 


10 


3 




5.92 


1420 


5 


0.5 


-0.13 


2.83 


2326 


5 


0.080 




2.75 



Table 2. Calibration errors cr ca i, zero level errors ctqi and baseline corrections from the original experiments at 408, 
142 and 2326 MHz; extragalactic background for the 408 and 1420 MHz maps have been estimated by Lawson et al. 
(1987), for the 2326 MHz map in this work following the same procedure. 



Fig. 7. Map of the differences Df3 between the spectral index evaluated with the destriped and the original maps. 
Features corresponding to stripes are evident. 



Fig. 8. Distribution of the differences (3, 



408,1420 



1420,2326 



The binsize is 0.01. 



for one data set, with no improvement. In the overlapping 
region, they set the spectral index to a weighted average 
of /3408,i420 and Ato8,2326; the choice of the weights ensures 
that the three different sky regions smoothly merge into 
each other. By convolving the data sets to 10° resolution, 
the stripes and discontinuities effects were averaged out. 
In this work we are interested in the high resolution infor- 
mation of the spectral index; after the destriping proce- 
dure discussed in Sect.|3J we investigate to what extent the 
discontinuities in the a and (3 maps could be accounted 
for within the quoted measurements systematic errors of 
the three surveys (see Table|5J), i.e., gain and zero-level er- 
rors. One should be aware that these uncertainties include 
a number of effects related to the details of the scanning 
strategy and calibration procedure adopted in the observa- 
tions. All three surveys are obtained by merging together 
different sky patches with zero levels somewhat different 
from one another. In the merging phase of the map pro- 
duction and as a-posteriori analysis, attempts were made 
to minimise discontinuities between different patches (e.g., 
Reich & Reich 1988). In addition, these surveys have been 
calibrated by convolution to sky horn absolute measure- 
ments, typically obtained at fixed declinations. With in- 
creasing distance from the absolute calibration data, the 
uncertainty in the zero level could increase due to resid- 
ual baseline gradients. Therefore using zero-level errors 
as "baseline errors" is equivalent to the assumption that 
the underlying constant component is of the same order 
as the total uncertainty. Although not ideal, we follow 
this assumption since it is the only practical possibility at 
hand. A deeper analysis would require the availability of 
the original raw data of each survey to perform a complete 
reanalysis, a task which is outside the scope of this work. 

We searched for the combination of six parameters 
(three pairs of parameters, calibration factor g ca \ and zero 
level correction 8q\, each for each frequency map) to rescale 
the maps to account for their zero level and calibration 
errors with the aim to minimize the difference across the 
discontinuities. Of course, the assumption that these three 
couples of parameters are constant over the entire survey 
is an over-simplification; in fact it is very hard, if not im- 



possible, to recover the real pattern of the variations of 
these parameters from the considered data products (see 
below for an attempt in this direction). The six parame- 
ters influence the (3 and a values through the optimized 
temperatures of the maps: 



T, 



v opt 



T v g cal 



301 • 



(9) 



We adopted the following x 2 -like form of the functional to 
be minimized: 



N 

E 



4 d « + 

(Q!2d(«) 



«3d0)) 2 



(<*2u(i) 



a 



fit 

a 3u 



(0 
i)) 2 



.(0 



(10) 



where u and d refers to the up and down discontinuities 
respectively, the low indices 2 and 3 refers to j3 and a eval- 
uated by exploiting two or three frequency maps, and the 
sum is carried out over i — th azimuthal position of all the 
N pixels of the top and bottom border lines of the overlap- 
ping region. The minimization has been carried out by us- 
ing the MINUIT package of the CERN library 6 . Different 
approaches have been tested in this minimization: we left 
all the parameters free, we kept fixed each of them (e.g. 
5q\ = and g ca \ — 1 or viceversa for a given frequency 
map) while leaving the others free; we varied each of them 
at one time keeping fixed the others, and finally, we var- 
ied the parameters c5oi and <?cai of only a given frequency 
map by keeping the "original" parameters at the other fre- 
quencies 7 . By observing that two frequency maps (at 408 
and 2326 MHz) show, before destriping, almost vertical 
striations due to the scanning direction, we implemented 
also a modified version of this minimization technique: 
instead of searching for six overall parameters over the 
whole maps, we searched for six azimuthal (e.g. as func- 
tion of the right ascension) profiles of the three pairs of <5oi 
and 5 ca i by implementing, at each i — th azimuthal posi- 
tion, an overall minimization of the discrepances between 

6 http://cern.web.cern.ch/CERN/ 

7 Of course the choice <? ca i = 1 and So 
the "original" map. 
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Fig. 9. Differences in f3 in the case of statistical error only (a) and systematics included (b); pixel number [0-1024] 
corresponds to right ascension [—180°, 180°]. The curves show the up (al and bl) and down (a2 and b2) discontinuities 
before (green) and after (purple) minimization. It is worth to remind that the values A/3 and Aa before minimization 
depend on the kind of errors included in the analysis, via their influence on the fit in the case we use data at all the 
three frequencies. 

Fig. 10. Differences in a in the case of statistical error only (a) and systematics included (b). The curves show the up 
(al and bl) and down (a2 and b2) discontinuities before (green) and after (purple) minimization. See comments for 

Fig. m 



statistical statistical + systematic 





errors only 


errors 


5^1,408 


0.877 


0.973 


^1,1420 


1.100 


1.066 


Seal, 2326 


0.827 


0.793 


^01,408 


3.432 K 


0.676 K 


^01,1420 


-0.0410 K 


-0.0267 K 


^01,2326 


0.0366 K 


0.0357 K 


X 2 /d.o.f. 


14.381 


0.204 



Table 3. Calibration factors and zero level corrections 
for the three maps as recovered after the minimization for 
statistical errors only and with systematic errors includes. 
The x 2 /d.o.f. is also shown. 



the values obtained for (3 and a by using three frequency 
maps or two combinations of only two frequency maps for 
all the pixels at the same i — th azimuthal position but at 
the different declinations where three frequency maps are 
available (~ 90 pixels at each i — th azimuthal position). 
This approach should in principle better account for pos- 
sible azimuthal variations <5oi and g ca \ over the maps in the 
case in which they do not show a significant declination 
dependence. Unfortunately, this method, although quite 
effective in reducing horizontal discontinuities, introduces 
undesired vertical discontinuities, suggesting for complex 
patterns of 8q\ and g ca \ whose recovery cannot clearly be 
derived by using such kinds of minimization approaches. 

We then adopted the overall minimization approach 
described by the eq. I|10|) searching for the optimized set 
of the above six parameters. 

Of course, the best result from the point of view of the 
reduced x 2 an d of the minimization of discontinuities, is 
obtained leaving the all six parameters free to vary. On 
the other hand, it is clear that we need to keep them in- 
side their corresponding experimental quoted errors. This 
requirement can be relaxed for the parameter g ca \ refer- 
ring to the map at 2326 MHz because of the uncertainty 
related to the polarization, as discussed in the previous 
subsection, which could account for variations of « 10%. 
Moreover, the 2326 MHz data have been calibrated using 
2 GHz data (Bersanelli et al. 1994) scaled with a spectral 
index of 2.75; the relaxed requirement can take into ac- 
count also the bias introduced in the determination of the 
spectral index distribution by the calibration procedure. 



We carried out the above minimization with <5oi and 
<7cai inside their quoted 2cr errors for the frequency maps 
at 408 and 1420 MHz, and with <5oi (ffcai) inside its quoted 
2cr error (5a error) for the frequency map at 2326 MHz. 

We have repeated this analysis by including or not 
the contribution of systematics in the error budget, as 
discussed in the previous subsection. 

In addition, we have tested that the results do not 
depend on the adopted range of allowed parameter val- 
ues when the systematic errors are included. On the other 
hand, neglecting the systematic errors, the minimization 
tends to give a best fit value of g ca \ at 1420 MHz out of 
its quoted 2a error. In this case, we have then repeated 
the minimization by constraining g ca \ at 1420 MHz within 
different ranges (within 1, 2 and 5 a) and keeping the al- 
lowed ranges for the other variables as before. We find an 
impact of a few percent on average and always less than 
10% (with a change less than 1% for half of the pixels) on 
the maps of /3 and a. In Table we show the six parame- 
ters after minimization. The results obtained by including 
or not systematic errors are in agreement within la error 
in each variable of the fit. The parabolic errors quoted for 
each recovered parameter from the minimization code are 
much smaller than the corresponding experimental uncer- 
tainties and discrepancies found including or not system- 
atics, and we avoid to report them, being clearly "unphys- 
ical" 8 . As expected, the x 2 /d.o.f. is larger (smaller) than 
unity when systematic errors are neglected (included) in 
the analysis. It is worth to point out that the parameters in 
Table 13 have to be considered not as the "real" correction 
to the experimental values of baseline and calibration of 
each survey, individually considered, but as the result of a 
procedure devoted to obtain more physically meaningful 
distributions of spectral index and normalization factor 
from the joint exploitation of the three surveys. This is 
particulary clear in the case of the 2326 MHz map where, 
to allow the comparison of a linearly polarized survey with 
the intensity data, the recovered calibration parameter has 
been found well out the experimental values. 

The differences, A/3 and Aa, in correspondence to the 
discontinuities, before and after having applied the de- 
scribed minimization neglecting or including the system- 
atics, are shown in Figs. l9landl 101 Note that the differences 
are evaluated between the two adjacent lines across each 

8 Clearly, in the maps of the final error on (3 and a the con- 
tribution of these parabolic errors is also negligible. 
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discontinuities, while the minimization is carried out with 
the values of evaluated with two or three frequencies 
in each border line of the overlapping region. The algo- 
rithm has the same qualitative effect if we consider only 
statistical or also systematic errors. 

4.3. Results 

The resulting maps of 0, a and their errors are shown in 
Figs-EHSl 

From Figs. I§1 and 1101 it is clear that discontinuities do 
not completely disappear; they are reduced, but the mini- 
mization procedure is not able to account for the full com- 
plexity of the real experimental drifts of the six parame- 
ters. The increasing of erg and a a with declination are 
explained considering the correction factor c a appearing 
in Eq. J5J) to account for non-equal area pixels. 

The distributions of the spectral index are shown 
in Figs. El an d El we evaluated the average spectral 
index and the dispersion before and after minimization, 
as summarized in Table the effect of minimization is 
to slightly change both parameters. The values well agree 
with those found by Giardino et al. (2002) of ,$480,1420 = 
2.78 ±0.17 and /3 48 o,2326 = 2.75 ±0.12 . The minimization 
procedure does not alter significantly the pixel by pixel 
values of 0: typically, a variation greater than 0.1 occurs 
for a few % of the pixels. 

We have estimated the contribution of point sources 
on the 8 distribution as follows: we subtracted the point 
sources from the destriped maps (with the same technique 
explained in Sect. 0J) and performed all the same steps on 
the analysis. The main results is a sensible reduction of 
the low spectral index tail in the distribution 9 ; on the 
other hand, when extended regions are considered, this 
procedure can sensibly smoothing the outer bordes of the 
structure. 

The incidence of free- free diffuse contribution has been 
estimated as follows. We scaled the WMAP free-free tem- 
plate (Bennett et al. 2003) at 408, 1420 and 2326 MHz 
and subtracted from the destriped maps. We obtained a 
variation of larger than 0.1 for less than 8% of the pix- 
els; this difference is inside the deviations illustrated in 
Table 0] except for a few percentage of the pixels. Thus, 
we can conclude that the diffuse free-free emission does 
not sensibly affect the distributions obtained in this work. 

We can say that the diffuse Galactic maps resulting 
from this work are representative for synchrotron spectral 
index distribution and synchrotron normalization factor. 
Point sources contribution is limited to a small number 
of pixels that could be isolated and removed to extrapo- 
late the maps at high frequency and obtain synchrotron 
emission templates. 



The high /3 tail comes from extended sources not removed 
by the source selection method. 



5. Discussion and conclusions 

We analysed the contribution of systematic uncertainties 
in the diffuse Galactic emission surveys at 408, 1420 and 
2326 MHz, cleaning the striping effect arising in combina- 
tion with the scanning strategy. 

We used the low frequency destriped maps to obtain 
the spectral index and normalization factor distributions 
over the whole observed sky. 

We performed a minimization technique over the cali- 
bration and zero level errors of the surveys to reduce the 
discontinuities in the a and maps in correspondence to 
the borders of the 1420 and 2326 MHz maps. The artifi- 
cial lines of discontinuities are reduced but not eliminated, 
suggesting for a complex pattern of Sqi and <? ca i. 

We discussed the spectral index distribution recovered; 
by estimating the contribution of point sources and diffuse 
free-free emission, we concluded that the spectral index 
distribution is representative for the diffuse Galactic syn- 
chrotron emission between 408 and 2326 MHz. 

The new results of the WMAP experiment give a quite 
precise picture of the Galactic emissions between ~ 23 
and 94 GHz IjBennett et al. 2003|> . The theory predicts a 
break in the synchrotron spectrum of the Milky Way at 

22 GHz (|Voelk 1989H . observed in the WMAP data. In 
the synchrotron spectral index map between 0.408 and 

23 GHz produced by the WMAP team, emission from 
the Galactic plane is dominant (the North Galactic Spur 
is much less evident than at radio frequencies), since, with 
increasing frequency, the flat components of the spectrum 
become more important than the steep components. For 
example, the fact that in the WMAP data the North 
Galactic Spur is much less evident than at radio frequen- 
cies shows clearly that any template of the Galactic emis- 
sion is frequency dependent. This puts more interest in 
the comparison between the new data and those derived 
from radio surveys. 
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statistical 


statistical + systematic 






errors only 


errors 


before minimization 


f3±op 


2.742 ±0.145 


2.734 ±0.124 




% pixels with j3 < 2 


0.031 


0.015 




% pixels with /3 > 3.5 


0.066 


0.05 


after minimization 




2.744 ±0.122 


2.695 ±0.120 




% pixels with j3 < 2 


0.022 


0.009 




% pixels with /3 > 3.5 


0.012 


0.012 



Table 4. Mean values and dispersion of the spectral index distributions. The percentage of pixels with /? < 2 and 
> 3.5 is shown. 
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Fig. 13. Maps of the synchrotron spectral index errors <jp after minimization. Top: the analysis has been carried out 
with temperature statistical errors only. Bottom: systematic errors have been included. 

Fig. 14. Maps of the normalization factor errors a a after minimization. Top: the analysis has been carried out with 
temperature statistical errors only. Bottom: systematic errors have been included. 

Fig. 15. Analysis with statistical errors only. The histogram of the spectral indices before (a) and after (b) minimization 
are shown together with their average and dispersion. 

Fig. 16. Analysis with statistical and systematic errors. The histogram of the spectral indices before (a) and after (b) 
minimization are shown together with their average and dispersion. 
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